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Abstract 

We provide an explanation of the spin-freezing transition recently observed 
by Chou et al. (Phys. Rev. Lett. 71, 2323 (1993)) in La 2 - x Sr x CuO A 
for x ^ .02. We propose that topological excitations of the 2D Heisenberg 
quantum antiferromagnet having non-coplanar transverse components have a 
pair-interaction energy that qualitatively and quantitatively agrees with the 
observed values of spin-freezing temperature as a function of doping. 
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I. INTRODUCTION: 



One of the avenues that may aid in the understanding of the properties of high- 
temperature superconductors involves examining the magnetic correlations that exist in the 
Cu02 planes when doped with carriers. Experiments may be employed to investigate the 
local magnetic character of these systems - some of the most useful are nuclear quadrupole 
resonance (NQR) studies where the nucleus provides a local probe |T],0], and from this in- 
formation many features of the magnetic ordering and/or fluctuations may be identified. 
Further, through resonance frequencies, one can measure the coupling between the nucleus 
and the electric field gradient surrounding it, thus providing insight into the electronic charge 
distribution of the system. 

The parent compound of the Bednorz-Miiller high-T c compound, viz. La^CuO^ is an 
antiferromagnetic (AFM) insulator which undergoes a Neel-ordering transition at 
~ 300 K. When the compound has trivalent La substituted for with divalent Sr for low 
doping levels, viz. less than 1% Sr, the AFM ordering persists with the Neel transition 
temperature (TV) rapidly decreasing with increasing Sr concentration. Recent l39 La NQR 
experiments || on lightly doped La2- x Sr x CuO^ (x < 0.02) have suggested that an inter- 
esting addition to the AFM structure takes place. To be specific, they have measured data 
consistent with the coexistence of the AFM phase with long-ranged order and a spin-glass 
phase, and have noted that this latter state can be interpreted to be composed of frozen 
transverse (i.e., perpendicular to the direction of the staggered magnetization) spin compo- 
nents. The state may be unambiguously assigned to be that of a spin glass via the stretched 
exponential recovery of the nuclear magnetization ||. A schematic diagram of the low x 
phase diagram is shown in Fig. 1, where a line separating the pure Neel state and the 
spin-glass coexisting state is denoted by Tf, we shall refer to this line as the spin-freezing 
temperature. The focus of this paper is to (i) suggest the spin texture of this coexisting 
state, and (ii) show that this spin texture predicts an onset of coexistence at Tf consistent 
with experiment. 
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To begin, let us summarize the experimental results. Chou et al. || have measured the 
139 La NQR for different dopings and temperatures in lightly doped La2- x Sr x CuO^ for low 
temperatures. They measured the 139 La NQR frequency spectrum, and the relaxation times 
for the system's magnetization to return to equilibrium. The splitting of the 2vq NQR 
peaks due to the internal field associated with the AFM ordering of the spins at the Neel 
temperature led to a doping-dependent ordering temperature T N (x) - this determination 
of Tn(x) is in agreement with other measurements, such as (i) their observations of the 
spin-lattice relaxation rate, and (ii) that found using static susceptibility measurements ||. 
Further, a large enhancement of this splitting was found at very low T, viz. at temperatures 
less than ~ 30 K, with the enhancement increasing with increasing x. Such behaviour could 
be caused by an alteration of the underlying spin texture at this temperature. A change 
in the spin state may also be inferred from measurements of the relaxation rate of the 3z/q 
transition - below ~ 70K there is a large enhancement of this rate. The rate continues to 
increase until it rapidly drops off at very low temperatures, viz. temperatures around 5 - 
16 K. Using this low temperature behaviour of the relaxation rate, Chou et al. || identified 
the maximum of the NQR rate with the spin-freezing temperature Tf. As a function of 
doping it was determined to be approximately given by Tf(x) ~ (815K)x, and it is this 
linear relation that we have sketched in Fig. 1. Above 30K (which may be interpreted to be 
the charge localization transition temperature |J) the nuclear magnetization relaxed with a 
multiexponential decay, while below this temperature a dramatic change was noted in that 
now the relaxation of the magnetization was found to be described by M(oo) — M(t) oc 
exp[— (i/T*) 1 / 2 ]. This stretched exponential decay has been found in many systems studied 
by NQR (see, e.g. Ref. ||), and can be associated with a spin-glass ordering. However, the 
staggered moment of the AFM phase does not disappear, as indicated, e.g., in (i) the zero- 
temperature extrapolation of the 2vq splitting to the undoped value, or (ii) the continuous 
change of the splitting as Tf is traversed (see Fig. 1(c) of [|] for the data corresponding 
to both of these results), and thus the spin glass ordering is found to coexist with the 
long-ranged AFM order. 
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Similar experimental results were found in both other NQR work at x = .012 and 
.015, as well as in electron paramagnetic resonance (EPR) studies of the same system at 
x = .009 by Rettori et al. |J (see their Fig. 12) - his latter study again found an increase in 
internal field observed in the —7/2 — > — 5/2 transition. 

Clearly, it would be desirable to explain the above-mentioned behaviour of the spin 
texture, and associated experimental measurements, and that is the focus of this paper. 
However, we wish to stress that this phenomenon is certainly not limited to the weakly 
doped Bednorz-Miiller high T c compound, nor to AFMs. For example, previous exper- 



imental work on metallic spin glasses ||, and related computational studies [|10j on the 
moderately frustrated three-dimensional Heisenberg model, aids in the understanding of the 
high T c system we are studying. The work on the moderately frustrated Heisenberg model 



leads to the phase diagram shown in Fig. 2 of Ref . |H| , and for dopant levels less than 25 
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is similar to that of our Fig. 1. In particular, the linear dependence of T xy is the same. 



The authors of Ref. [T(J identified the coexisting freezing of the transverse spin degrees of 
freedom (transverse in the sense of being perpendicular to the already present ferromagnetic 
moment) and long-ranged ferromagnetic order at the spin-freezing temperature T xy ; this 
latter ordering is not associated with a phase transition, only involving a short-ranged freez- 
ing of spins, and does not eliminate the ferromagnetic ordering. Unfortunately, these authors 
did not attempt to come to an analytic understanding of this similar linear dependence of 
the spin-freezing temperature on doping. Nor did they attempt to justify the magnitude 
of the slope of this pseudo-phase boundary. That is precisely what we will present in this 
paper, albeit for at a doped two-dimensional AFM. 

Our paper is organized as follows. In §11 we summarize our previous work |TT|JT2| on the 
spin texture of singly and weakly doped CuO<i planes; this involves studies of a non-coplanar 
spin texture where each Sr 2+ dopant produces a spin state reminiscent of a singly charged 
skyrmion. This section also serves to introduce the effective Hamiltonian that we will use to 
describe these doped planes. In §111 we present an analytic calculation of the two skyrmion 
interaction energy, and then use this to suggest the dependence of Tf(x) on doping x. In §IV 



we describe numerical work examining the same two skyrmion interaction energy problem, 
and which reaches the same qualitative conclusion as in §111, viz. an interaction energy oc x, 
as well as providing a quantitative estimate of the slope. Then, we conclude this section by 
extrapolating the above work to determine the dependence of the spin-freezing temperature 
on doping via both a comparison to Chou et al. [f5|, and a simple counting of the degrees of 
freedom of the two defect system. Finally, in §V we summarize how this and other work on 
the weakly doped planes leads to a clear description of the novel spin texture of this system. 



II. SPIN TEXTURE OF WEAKLY DOPED PLANES AT LOW T: 

Our model of the doped state of a CuOi plane at low T is based on a collection of partially 
delocalised holes. It is assumed that at low temperatures the carriers (viz. O holes) are not 
mobile, but rather are confined to some small region of a plane - this is entirely consistent 



with the localisation transition seen in resistivity measurements between 50 and 100 K [JT3 
as well as with the localisation transition inferred from the NQR measurements of Chou || . 
Physically, this localisation can be considered to arise from the electrostatic attraction of 
the holes to a divalent Sr 2+ ion. 

A study of the magnetic spin texture of the ground state of the system for one hole was 



previously discussed by one of us [11], within the framework of the t — J model. There, 



the hole was only free to move around one square plaquette, and it was found that the 



distorted AFM spins' state could be described, in a semi-classical sense fl4}[ , via a non- 
coplanar spin state with a rotational twist generated by the motional current of the hole 
circulating around the plaquette. Since the hole can circulate in either the clockwise or 
counter-clockwise directions, the ground state is two-fold degenerate. The combination of 
the above results allows one to infer that the spin texture is analogous to a singly charged 
skyrmion found in classical chiral field theory |I5]] . A numerical study JIB] found similar 
results in that the symmetry of the ground state in their work is identical to that of the 
skyrmion. 
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In order to model many such partially localised carriers an effective Hamiltonian based 
on only the spin degrees of freedom was recently proposed by two of us JTIJ . Again, based 
on the semiclassical representation of the spins, a Hamiltonian thought to mimick the effect 
of a hole confined to a single plaquette was used: 

H= JES,-S,-^{(S 1 -S 2 xS 3 ) 2 

+ (S 2 • S 3 x S 4 ) 2 + (S 3 • S 4 x Sx) 2 + (S 4 • Si x S 2 ) 2 } . (2.1) 

The first term is the usual Heisenberg interaction representing the AFM background. J is 
the exchange coupling and (ij) denotes nearest neighbours in the plane; all near-neighbour 
pairs are summed over. The second term is introduced to generate a spin distortion that 
is the same as that produced by a hole circulating around the plaquette, the corners of the 
plaquette being denoted by 1, 2, 3 and 4. It was motivated by the work of Wen, Wilczek, 
and Zee O] on chiral spin liquids, and we shall refer to this term as the chiral-plaquette 
interaction. The prefactor of the interaction, D/S A , includes the l/S 4 factor to ensure that 



the Heisenberg and one-plaquette interactions both scale as S 2 . In Ref. [12] it was found 
that the interaction strength D must exceed ~ 2.2 J to be able to induce a singly-charged 
skyrmion ground state. 

In order to model a plane doped with many O holes, one simply sums the chiral-plaquette 
interactions over all plaquettes having a Sr 2+ ion above them; we distribute these defects 
randomly in the plane, and the resulting system is necessarily spatially inhomogeneous. The 
results of Monte Carlo studies of the spin-spin correlation function [18| for large lattices at 



temperatures greater than J/2 were consistent with the empirical relation 

r 1 (x 1 T) = r 1 (x,o)+r 1 (o,T) (2.2) 

for the spin correlation length £(x,T) found experimentally by Keimer et al. [0. Note that 
in order for Eq. (2.2) to be obtained, it was necessary to use D ~ 3 J. Further, in Ref. 
T2^ we showed that a mean field model of this effective Hamiltonian was also consistent 



with Raman scattering experiments ||20|| . Now, in this paper, we extend our study of this 



effective Hamiltonian and demonstrate that the transverse spin-freezing temperature Tf 
may be quantitatively determined. 



III. SEMICLASSICAL EVALUATION OF Tf. 

Here we present a discussion that leads to the analytic prediction for the dependence 
of the spin-freezing temperature on doping. The theory is based on a calculation of the 
interaction energy between two skyrmions, and is performed utilizing the semiclassical ap- 
proach of Shraiman and Siggia . Two of us |2l| have previously used a similar approach 
to calculate the interaction energy between two spin polarons generated in the frustrated 



bond model ||22|.|23| , a calculation that was also verified via exact diagonalization numerical 



work |gT|. 



To begin, we introduce the spin texture associated with one skyrmion, where the spins are 
described by a continuum field which is chiral. The direction of the staggered magnetization 
is denoted by a chiral vector field f2 satisfying 

£i £1 = 1 (3.1) 

at every point in the plane. For a singly charged skyrmion (Q=± 1) placed at the origin of 
the coordinate system in use, and for one choice of phase (see below), where the staggered 
magnetization is chosen to be in the positive z direction at infinity, the chiral field is found 
to be HI 



A 2Xx 2Xy A r 2 -A 2 A , . 

X z + r z X z + r 2 r A + X 2 

where r = \Jx 2 + y 2 is the distance in the two-dimensional plane from the centre of the 
skyrmion, and A is the size of the skyrmion. That A is the size of the skyrmion is seen 
from the z component, viz. for r < A the direction of the staggered magnetization is 
reversed; formally, A sets the length scale in the continuum problem. Further, upon changing 
the topological charge Q from 1 to - 1, the handedness of the twisting of the staggered 
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magnetization is found to change. This is seen in the x and y components: note that 
Cl x ± iCly oc x ± iy, and thus the twisting of the staggered magnetization is described 
by a vector which circulates in the counter-clockwise (Q = 1) or clockwise (Q = - 1) 
directions. Lastly, note that these two components of the spin field represent the transverse 
degrees of freedom and in this paper we shall suggest that it is these components that lock 
together, due to the mutual interactions between skyrmions mediated by the Heisenberg 
AFM superexchange, leading to the spin freezing observed experimentally. A schematic 
representation of the spin texture for the Q = 1 skyrmion is shown in Fig. 2. 

We place one skyrmion at r — 0, and then a second at r = (£,0). A charge of Q = 
1 is chosen for the first, and the charge of the second will be determined by minimizing 
the interaction energy. The latter may be obtained from extrapolating the semiclassical 
theory ||14| to describe this system. There, in the classical approximation, Shraiman and 



Siggia derived that the interaction between a dipolar spin distortion P a , such as that which 
would be produced by a ferromagnetic bond in an AFM background (here a = x or y labels a 
direction in the plane, such as the direction of the ferromagnetic bond), and the surrounding 
spins is given by 

Edassical = ~92 ^ P « " X d afr) ■ (3-3) 



The interaction strength g 2 is a phenomenological constant [|l4j] which in the t ^> J limit 
is of O(J). The ground state for the spin texture is then found by minimizing the sum of 
this energy and the usual Heisenberg AFM interactions between spins, the latter being most 
conveniently represented by the nonlinear sigma model. We cannot associate a dipole mo- 
ment with the spin texture of a skyrmion, although it can be represented as a superposition 
of dipolar backflow spin distortions induced by mobile holes subsequently localised onto a 
single plaquette |M| . A more direct way to represent the interaction of the background spin 
texture with the distortion induced by the skyrmion is to replace 

P a - Cl x d a Cl (3.4) 



S 



where the staggered magnetization in this equation is that of the skyrmion. This is the 
analogue of the definition of the semiclassical dipole moment given by Shraiman and Siggia 
p5| for mobile holes. Then, we posit that one can represent the interaction between two 
skyrmions as 

E mt = -g 2 £ (A« x d a &W) ■ (tt 2) x d a Cl^) (3.5) 

a=x,y 

where the superscripts 1 and 2 denote each of the two skyrmions. (Here we shall assume 
that each of the spin fields describing a skyrmion are those produced by systems with only 
one defect. Thus, no back-reaction effects are included - these will be considered in the 
next section.) As an estimate of this energy we evaluate E int at either of the points in space 
where either of the two skyrmions are located. In the limit of large £ we obtain 

E int = -8J/f, £ -» oo (3.6) 

when the second skyrmion is chosen to be a Q = - 1 spin distortion, as this leads to the 
minimum energy. We note that the second skyrmion can be placed anywhere in the lattice 
with the same result (viz. Eq. (3.6)) being obtained as long as the second skyrmion has Q 
= - 1. However, an arbitrary phase of the skyrmion must be introduced (one that is not 
included in Eq. (3.2) - see, e.g., Ref. |1| for such a spin field), and then fixed via minimizing 
E int with respect to this angle. (The interaction energy between two skyrmions with the 
same topological charge scales as l/£ 4 in the large £ limit. As we will see below, this would 
lead to Tf oc x 2 , which indeed is not observed experimentally.) 

For holes randomly placed on a two-dimensional lattice one expects that the doping x 
would be given by x = l/£ 2 . This is consistent with low temperature neutron scattering 
studies of Keimer et al. [|19j where the correlation length £(x,T) was found to be behave 
approximately as £(x, 0) ~ l/\/x. Thus, expressing the two-skyrmion interaction energy as 
a function of doping, we finally have 

E int = -8Jx. (3.7) 
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Now, we suggest that one may associate the spin-freezing temperature Tf with this 
interaction energy; to be specific, ksTf oc |-Ei„t|. This is consistent with the notion that 
there is only short-ranged ordering of the transverse spins (similar to the results of Ref. 
||10|| ) and the temperature at which the freezing occurs is such that at higher temperatures 
thermal fluctuations would destabilize a spin texture having these transverse degrees of 
freedom fixed. Of course, the energy associated with when the transverse degrees of freedom 
would become frozen into their ground state involves many skyrmions, and thus we can 
only consider our treatment to be valid in the low doping limit. Noting that the spins are 
S — 1/2 quantum spins, and thus we must include the renormalization factor Z c accounting 
for quantum fluctuations, we finally have that 

k B T f oc 8JZ c x (3.8) 



where Z c ~ 1.2 pq| . This relationship demonstrates that if the transverse degrees of freedom 
which freeze are associated with the many skyrmion spin texture produced by randomly dis- 
tributed defects, a linear dependence of Tf on doping is expected. This is entirely consistent 
with the experiments of Chou et al. [^J. 



IV. NUMERICAL EVALUATION OF Tf. 

We have just seen that the freezing of the transverse degrees of freedom of two skyrmions 
separated by a distance 1/y^ would lead to a spin-freezing temperature linear in doping 
x. In this section we present a numerical procedure to obtain the quantitative dependence 
of Tf on the doping again in the low temperature and lightly doped regime. This work 
is required since the relationship found in the previous section, viz. Eq. (3.8), did not 
allow for (i) the back-reaction, and subsequent alteration, of a skyrmion's core, and thus 
must be an overestimate of the interaction energy, and (ii) relies on an assumption for #2, 
the interaction strength. Here we wish to eliminate these limitations and determine a very 
accurate numerical relationship for Tf on x. 
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To establish a relationship between Ei nt and x we will consider a lattice of classical spins 
with two skyrmions placed somewhere in the lattice a distance I lattice spacings apart. 
We calculate the interaction energy for such configurations by evaluating the ground state 
using Eq. (2.1) as the Hamiltonian (we have done the calculations for S — 1, and properly 
convert to S = ~ below). In Ref. [|r2] it was found that the constant D was required to be 



2.5 < D/J < 3.3 in order to explain both the Raman scattering data j20[ and the neutron 
scattering experiments consistent with Eq. (2.2). Here we perform the evaluation of Ei nt 
for D/J = 2.5,3, and 3.3 in Eq. (2.1) to see how Ei nt varies with D/J. 

The method employed in the evaluation of Ei nt is straightforward and involves studying 
LxL lattices with open boundary conditions with two chiral-plaquette interactions included 
on two plaquettes, and then extrapolating to L — > oo; for concreteness, we explain the 
method by considering Fig. 3. In this figure we have shown a 6 x 6 lattice and have labeled 
two plaquettes Si and S2, where the separation between the centres of these plaquettes 
correspond to I = 2. Firstly, we find the single skyrmion energy for finite L by eliminating 
the chiral-plaquette interaction on S2, and then find the ground state spin configuration 
and associated energy (with open boundary conditions) - denote this energy by Ei(L,£), 
where the inclusion of £ in E\ indicates that the skyrmion is off centre on the lattice. Then, 
the chiral-plaquette interaction on S2 is reintroduced and the two-skyrmion energy may be 
found once the two-skyrmion ground state energy, denoted by E 2 (L,£), is known (again 
evaluated for open boundary conditions). All zero-force spin configurations were evaluated 



using a conjugate-force method similar to that employed in Ref. ||21|| . 

Consistent with the above semiclassical work, we always find that the ground state for 
two defects corresponds to opposite topological charges on the plaquettes. Then, in analogy 
to a two-particle binding energy, one defines 

E int (LJ) = E 2 (L,£)-2(E 1 (L,e)-E (L))-E (L). (4.1) 

where E (L) is the energy of a pure Heisenberg AFM, with open boundary conditions, for 
an L x L lattice. Finally, the large L limit is taken to evaluate the interaction energy. 
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An example of this extrapolation is shown in Fig. 4 for D/J = 3 and £ = S\/2, the latter 
corresponding to x — .0078. 

In Fig. 5 we have plotted \Ei nt \/J vs x, utilizing i = l/y/x, for all three ratios of 
D/J. We obtain 5.14a; ± 0.08, 5.00a; ± .08, and 4.74a; ± 0.08, for D/J = 3.3, 3.0 and 2.5, 
respectively. The linear relation found in Eq. (3.8) is also found in the numerical work, 
and thus validates the hydrodynamic description of the spin degrees of freedom employed 
in the previous section. Further, it is clear that Eq. (3.7) represents the D — > oo limit 
(viz. no alteration of the skyrmions' cores), and thus the energies found in this section are 
appropriately less than that of Eq. (3.7). It is encouraging that these energies do not change 
that much with D/J. 

To determine the quantitative behaviour of Tf(x) we must (i) transform the above ener- 
gies to represent those of S = | quantum spins, and (ii) determine the ratio of \E int \ to ksTf. 
The first point may be incorporated by multiplying the interaction energies by S(S + 1) Z c . 
The second is a nontrivial exercise, and to resolve it one can appeal to experiment. Chou 
et al. || investigated how best to describe the low-temperature maxima they observed in the 
NQR rate, and found that if one attempted to describe the maxima by a non-equilibrium 
freezing of the spins, with a correlation time r associated with an Arrhenius law, viz. 

r oc exp(E/k B T), (4.2) 

then the relaxation rate peaked at Tf with a related exponent 

(E/k B T f ) = 8.9 ± 0.5 (4.3) 

(see Fig. 3 of Ref. j3|). Clearly, the energy E in the exponent in Eq. (4.3) must be 
associated with an interaction energy, and thus we incorporate \E int \ into this equation. 
A theoretical explanation for this relationship follows from application of the equipartition 
theorem. Apart from continuous symmetries, there are three differing low-energy states of 
the single defect problem: Q = ± 1, as well as an excitation with zero topological charge 



(see Ref. [|TTJ for a discussion of this important low-energy excitation). Thus, there are 3 2 
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static degrees of freedom for the two defect system. However, since the spins are quantum 
spins, one should include a kinetic or dynamical component for the quantum fluctuations 
associated with each of these degrees of freedom. Then, counting \ksT for each such degree 
of freedom one expects 

Eint = 2 3 2 h B T f . (4.4) 

This is consistent with the experimental determination of Chou et al. given in Eq. (4.3). 
Then, we immediately have estimates for the spin-freezing temperature: 

T f ~ (805 K) x, for D = 3.3 (4.5) 
T f ~ (784 K) x, for D = 3.0 (4.6) 

Tf ~ (743 K) x, for D = 2.5 (4.7) 

where we have used J ~ 1550 K for the exchange constant. We anticipate error bars of at 
least 5% on the above values, the major uncertainties coming from Z c , and the exchange 
constant J. 

The experimental result, viz. Tf ~ (815 K) x, shows that for D being close to 3 
this two-skyrmion interaction model provides a very credible estimate of the spin-freezing 
temperature. We stress that the temperature and doping dependence of the spin correlation 



length was found to agree with Ref . [19| , and the Raman scattering spectra of Ref . [ 20| , when 
the interaction strength D ~ 3 was used. It is gratifying that quantitative aspects of all three 
experiments seem to be consistent with the spin texture found when the low-temperature 
doped planes are described by Eq. (2.1). 

V. DISCUSSION: 



We have considered a simple model of the ground state spin texture produced by a 
Sr 2+ defect and a single hole localised in its region [11|, and then utilized a model that 
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allows for the description of many such defects [0, viz. Eq. (2.1), where the chiral 



plaquette interaction is summed over all plaquettes having a divalent Sr above them. We 
have used a semiclassical theory, and a numerical simulation technique, to determine the 
interaction energy between two such skyrmion states. We find that the ground state of 
such a system corresponds to two skyrmions having opposite topological charges, and that 
then the interaction energy, analogous to a two-particle binding energy, is proportional to 
l/£ 2 , as in Eq. (3.6), where £ is the separation between defects. Then, we associate the 
doping x with this average distance via i = 1/y/x, consistent with the low-temperature 
spin correlation length determined by neutron scattering studies. The interaction energy is 
then found to scale linearly with doping, and if we associate the spin-freezing temperature 
with this interaction energy we can reproduce the same linear dependence that is observed 
experimentally. We have quantified this relationship by using the experimental value for the 
Arrhenius exponent given in Eq. (4.3). We can also present theoretical arguments that lead 
to values for the ratio of E to kgTf that give the same numerical value as in Eq. (4.3) (e.g., 
a simple equipartition theorem argument also gives a ratio of 9). 

Alternatively, the frustrated bond model of Aharony et al. |22| predicts an interaction 



energy which scales linearly as x (see Ref. PT|), but the spin texture of the ground state 
depends on the location of the frustrated bonds in the plane (e.g., see Ref. plR ). Also, 
the ratio of Eq. (4.4) would be 4 for the frustrated bond model, and not ~ 9, as found 
experimentally. However, it is not clear whether or not the experimental results of Chou 
et al. f3[ depend on probe frequency, and thus it is not clear if this particular kind of 
experiment could distinguish between these two models of the spin texture of the weakly 
doped planes. Since the important physics of both these models is that defects produce long- 
ranged spin distortions, and there is some evidence that such distortions are also produced 
by mobile holes at moderate doping levels |fT^ , p7| , it is encouraging that both models rely 
on spin deviations coupling the perturbing effects of the holes to the long-wavelength spin 
waves, thus further justifying our hydrodynamic approach. 

We have showed that for our model, with the interaction strength D/J ~ 3 in Eq. (2.1), 
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the dependence of Tf on x is reproduced quite accurately. Further, since our model with 
this same ratio of D/J also reproduces the neutron scattering measurements of the spin 
correlation length |T9| , |28| and the Raman scattering spectra |2(J, it seems likely that these 
novel topological excitations of the S = | quantum AFM do indeed play a role in forming 
the spin textures of the low-temperature, inhomogeneously doped Cu02 planes. 

We wish to thank Fangcheng Chou, Ferdinand Borsa, and David Johnston for valuable 
discussions. This work was supported by the NSERC of Canada. 
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FIGURES 



1. Schematic phase diagram of the weakly doped system, showing that at temperatures 
below Tf(x) a coexisting quantum AFM with long-ranged order and a spin-glass phase is 
observed. 

2. Spin texture of a skyrmion state, such as that given in Eq. (3.2), shown on a small 6x6 
cluster (here, Q=l, and A = .5). The spin state shown in the figure is found by taking Eq 
(3.2) and rotating all spins (in spin space) by n/2 along the +z direction. The spins at 
infinity are chosen to all point out of the page, and we have only displayed the projection 
of the spins onto the plane. Also, for ease of presentation, all spins on B sublattice sites of 
the background Neel state have been inverted (Si — > —Si), making the D = state be a 
completely polarized ferromagnetic state. The solid dots denote the lattice sites. 

3. Schematic diagram of a small cluster denoting two Sr defects above plaquettes Si and 
S 2 - 

4. Interaction energy defined in Eq. (4.1), in units of J, for D = 3 and a separation of 
defects of 8a/2 (viz., eight diagonal spacings), as a function of system size L for L = 20, 24, 
and 28, showing the extrapolation to the bulk limit. 

5. The magnitude of the interaction energy \E int \, in units of J, as a function of doping x 
for differing ratios of D / J. For each D/J, the data extrapolates linearly back to the origin. 
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